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Abstract 



^ . The color evaporation model simply states that charmonium production is 

1/^ I described by the same dynamics as DD production, i.e., by the formation 

I of a colored cc pair. Its color happens to be bleached by soft final-state 

ly-^ ' interactions. We show that the model gives a complete picture of charmonium 

O ■ production including low-energy production by proton, photon and antiproton 

. beams, and high-energy production at the Tevatron and HERA. Our analysis 

i includes the first next-to-leading-order calculation in the color evaporation 
model. 
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I. INTRODUCTION 



In a recent paper [Q] we pointed out that an unorthodox prescription for the production of 
rapidity gaps in deep inelastic scattering, proposed by Buchmiiller and Hebecker p[ , suggests 
a description of the production of heavy quark bound states which is in agreement with the 
data. The prescription represents a reincarnation of the "duahty" or "color evaporation" 
model. It is very important to study the validity of this approach as it questions the 
conventional treatment of the color quantum number in perturbative QCD. 

The conventional treatment of color, i.e., the color singlet model, has run into serious 
problems describing the data on the production of charmonium While attempts to rectify 
the situation exist ||^, the color singlet model has remained the standard by which other 
approaches are measured. Here we show that the color evaporation approach, which actually 
predates the color singlet approach, describes the available data well. After a discussion of 
the model itself, we show how the ratio of production by antiprotons to production by 
protons forms a quantitative test of the model. The model passes the test extraordinarily 
well. We then show how further quantitative tests can be performed using the the data on 
photoproduction and hadroproduction total cross sections. We finally show that the model 
successfully describes pr distributions at the Tevatron. 

II. THE MODEL 

Color evaporation represents a fundamental departure from the way color singlet states 
are treated in perturbation theory. In fact, color is "ignored". Rather than explicitly impos- 
ing that the system is in a color singlet state in the short-distance perturbative diagrams, 
the appearance of color singlet asymptotic states depends solely on the outcome of large- 
distance fluctuations of quarks and gluons. These large-distance fluctuations are probably 
complex enough for the occupation of different color states to respect statistical counting. 
In other words, color is a nonperturbative phenomenon. In Fig. |I| we show typical diagrams 
for the production of ^/'-particles using the competing treatments of the color quantum num- 
ber. In the diagram of Fig. |I|a, the color singlet approach, the ip is produced in gluon-gluon 
interactions in association with a final state gluon which is required by color conservation. 
This diagram is related by crossing to the hadronic decay ip ^ ^ gluons. 

In the color evaporation approach, the color singlet property of the ip is initially ignored. 
For instance, the ip can be produced to leading order by gg-annihilation into cc, which is 
the color-equivalent of the Drell-Yan process. This diagram is calculated perturbatively; 
its dynamics are dictated by short-distance interactions of range Ax ~ m^^. It does not 
seem logical to enforce the color singlet property of the ip at short distances, given that 
there is an infinite time for soft gluons to readjust the color of the cc pair before it appears 
as an asymptotic or, alternatively, DD state. It is indeed hard to imagine that a color 
singlet state formed at a range m^^, automatically survives to form a ip. This formalism 
was, in fact, proposed almost twenty years ago and subsequently abandoned for no 

good reason. 

The color evaporation approach to color leads to a similar description of bound and open 
charm production. In the color evaporation model 
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FIG. 1. Typical diagrams for (a) color singlet ip production and (b) color evaporation ip pro- 
duction. 
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where the cross section for producing heavy quarks, aec, is computed perturbatively. Dia- 
grams are included order-by-order irrespective of the color of the cc pair. The coefficients | 
and I represent the statistical probabilities that the 3x3 charm pair is asymptotically in a 
singlet or octet state. 

The model also predicts that the sum of the cross sections of all onium states is given 
by Eq. ([^). This relation is, unfortunately, difficult to test experimentally, since it requires 
measuring cross sections for all of the bound states at a given energy. 

Other approaches similar in spirit can be found in Refs. @] and . The color evaporation 
approach differs from Ref. the formalism of Bodwin, Braaten and Lepage (BBL), in the 
way that the cc pair exchanges color with the underlying event. In the BBL formalism, 
multiple gluon interactions with the cc pair are suppressed by powers of v, the relative 
velocity of the heavy quarks within the tp. The color evaporation model assumes that these 
low-energy interactions can take place through multiple soft-gluon interactions; this implies 
a statistical treatment of color. 

The color evaporation model assumes a factorization of the production of the cc pair, 
which is perturbative and process dependent, and the materialization of this pair into a 
charmonium state by a mechanism that is nonperturbative and process independent. This 
assumption is reasonable since the characteristic time scales of the two processes are very 
different: the time scale for the production of the pair is the inverse of the heavy quark mass, 
while the formation of the bound state is associated to the time scale 1/Aqcd- Therefore, 
comparison with the ip data requires knowledge of the fraction of produced onium states 
that materialize as ip^s, i.e., 



0"^ — Pv'^o 



(3) 
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where is assumed to be a constant. This assumption is in agreement with the low energy 
data ||TO|JTl|| . We demonstrated in Ref. ^ that simple statistical counting estimates of p^, 
p-^, etc., accommodate all charmonia data to better than a factor of 2. 



III. THE TESTS 

We discussed several qualitative tests of the color evaporation picture in Ref. One 
such test that we did not mention is the polarization of produced charmonium. In the frame- 
work of the color evaporation model the multiple soft gluon exchange destroys the initial 
polarization of the heavy quark pair |jl2|. This fact is in agreement with the measurements 



of the ip polarization made in fixed-target vr N and pN ||TJ] reactions. The color singlet 
model fails to describe this feature of charmonium production since it predicts that ip^s 



are produced transversely polarized. The predictions in the literature for polarization in the 
BBL formalism have been somewhat controversial []TB|. 



One of the most striking features of color evaporation is that the production of char- 
monium at low energies is dominated by the conversion of a colored gluon into a. ip, as in 
Fig. |l]b. In the conventional treatment, where color singlet states are formed at the pertur- 
bative level, 3 gluons (or 2 gluons and a photon) are required to produce a ip. The result is 
that in the color evaporation model ip^s are hadroproduced not only by gluon-gluon initial 
states, but also via quark- ant iquark fusion. In the color singlet approach such diagrams only 
appear at higher orders of perturbation theory and their contribution is small. 

We can distinguish between the two pictures experimentally by comparing the production 
of charmonium by proton beams with that from antiproton beams. The color evaporation 
model predicts an enhanced ip cross section for antiproton beams, while the color singlet 
model predicts the same cross section for the production of -i/^'s, whether we use proton or 
antiproton beams. The prediction of an enhanced p yield compared to p yield at low energies 
is obviously correct: the ratio of antiproton and proton production of ^'s exceeds a factor 
5 close to threshold. (See Fig. Q.) In fact, it has been known for some time that ip^s are 
predominantly produced by qq states . Nonetheless, we should note that for sufficiently 
high energies, gluon initial states will eventually dominate because they represent the bulk of 
soft partons. This can be seen in Fig. ^ where the ratio gets close to unity for center-of-mass 
energies as low as 25 GeV. 

The merits of the color evaporation approach can be first appreciated by studying the 
data in calculation-independent ways. The factorization of the production of the cc pairs 
and the formation of onium states implies that the energy dependence and kinematic dis- 
tributions of the measured cross section for the different bound states should be equivalent. 
Moreover, in the approximation ttlc ^ mo this equivalence extends to the production of 
open DD pairs. In Fig. ^ we display charm photoproduction data for both open charm and 
bound state production with common normalization in order to show their identical energy 
behavior. In Fig. ^ we display charm hadroproduction data in a similar fashion. 

Further quantitative tests of color evaporation are made possible by the fact that the fac- 
tor is the same in hadro- and photoproduction. Once p^ has been empirically determined 
for one initial state, the cross section is predicted without free parameters for the other. We 
show next that color evaporation passes this test, quantitatively accommodating all mea- 
surements, including the high energy Tevatron and HERA data, which have represented a 
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FIG. 2. Ratio of the cross sections for the production of J/ip by proton and antiproton beams 
in the color evaporation model (solid line) and the color singlet model (dashed line) as a function 
of the center-of-mass energy. Data taken from Ref. [p^]. 
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FIG. 3. Photoproduction data |T|,|T9| clS Si function of the photon energy in the hadron rest 
frame, W^. The normahzation has been adjusted to show the similar shapes of the data. 



considerable challenge for the color singlet model. 

In Fig. ||we compare the photoproduction data with theory, using the NLO perturbative 
QCD calculation of charm pair production from Ref. |2^. The solid line is the NLO pre- 



diction for open charm production using the GRV HO distribution function and the scale 
/i = rric = 1.45 GeV. The dashed line is the prediction using the MRS A distribution function 
with the scale fi = 2mc = 2.86 GeV. We obtained these charm quark masses from the best 
fit to the data for each structure function. For the J/ip production data we employed the 
parameters used for describing open charm, and determined the fragmentation factor to 
be 0.50 using GRV HO, or 0.43 using MRS A. Note that the factor possesses a substantial 
theoretical uncertainty due to the choice of scales and parton distribution functions. We 
conclude the photoproduction of J/ip and DD is well described by the color evaporation 
model. This reaction determines the only free parameter, p^ 0.5. 

At this point the predictions of the color evaporation model for hadroproduction of ijj 
are completely determined, up to higher order QCD corrections. These can be determined 
by fitting the hadroproduction cross section of DD pairs with a global K factor. This 
factor is subsequently used to correct the ip prediction. In Fig. |^ we compare the color 
evaporation model predictions with the data. In order to fit the DD cross section with 
the NLO result for the production of cc pairs we introduced a factor K = 1.27 (1.71) for 
the GRV HO (MRS A) distribution function with the same scale p and charm quark mass 
determined by photoproduction. Inserting this K factor and the constant p^ obtained from 
photoproduction into Eq. (||), we predict ip hadroproduction without any free parameters. 
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FIG. 4. Hadroproduction data as a function of the center-of-mass energy, Ecm- Data taken 
from Refs. [21,^]. Once again, the normahzation has been adjusted to show the similar shapes of 
the data. 



7 



lot 



10' 



10^ 



1 — I — I I I 1 1 II 1 — I — I I I 1 1 II 1 — I — I I I 1 1 1 1 1 — I — I I I I Itl 



M =1.45; |Li = M ; GRV 

c c 

M =1.43; |Li = 2M ; MRS-A 




J I 



(GeV) 



FIG. 5. Photoproduction data |18| , |l9[ and the predictions of the color evaporation model 
at next-to-leading order as a function of the photon energy in the hadron rest frame, W^. The 
normalizations in this figure are absolute. 
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FIG. 6. Hadroproduction data |21,22| and the predictions of the color evaporation model at 
next-to-leading order as a function of the center-of-mass energy, Ecm- The curve for bound state 
production is an absolutely normalized, parameter-free prediction of the color evaporation model. 

We conclude from Fig. | that the color evaporation model describes the hadroproduction of 
ip very accurately. 

As a final test, we turn to pt distributions at the Tevatron. The CDF collaboration has 



accumulated large samples of data on the production of prompt ip, Xcj, and ip' [^. This data 
set allows a detailed study of the pr distribution of the produced charmonium states. Since 
all the charmonium states share the same production dynamics in the color evaporation 
model, their px distributions should be the same, up to a multiplicative constant. This 
prediction is confirmed by the CDF data, as we can see in Fig. |^. 

In order to obtain the theoretical prediction, we computed the processes g + g [cc] +g, 



q + q ^ [cc] + g, and g + q ^ [cc] + g at tree level using the package MADGRAPH p4 
This calculation should give a reliable estimate of the pt distribution at large values of pt- 
We require that the cc pair satisfy the invariant mass constraint on Eq. ([^). Our results 
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FIG. 7. Data from the CDF Collaboration |2S], shown with arbitrary normalization. The 
curves are the predictions of the color evaporation model at tree level, also shown with arbitrary 
normalization. The normalization is correctly predicted within a K-factor of 2.2. 

are shown in Fig. which required a K factor of 2.2 in order to fit the total ip production. 
As we can see from this figure, the color evaporation model describes the general features 
of the Pt distribution of the different charmonium states. Furthermore, we should keep in 
mind that the factor K is, in general, pt dependent and that it is usually larger at low pt 
Higher order corrections such as soft-gluon resummation are expected to tilt our lowest 
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order prediction, bringing it to a closer agreement with the data pH . 



IV. CONCLUSIONS 



The color evaporation model died an untimely death. In its current reincarnation its 
qualitative validity can be proven directly from the available experimental data by taking 
ratios of production by different beams and production of different particles. Moreover, we 
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have shown that its vahdity extends to the quantitative regime once we use next-to-leading- 
order QCD calculations for the photo- and hadroproduction of ip and DD pairs. 

The color evaporation model explains all available data on px distribution and energy- 
dependence of the cross section for the production of t/j in all the available energy range. 
Moreover, it sheds some light in the relation between the production of charmonium states 
and DD pairs. 
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